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Photoionization cross sections and partial ion yields of Xe and XeF 2 from Xe 3d5/2, Xe 3d3/2, and F Is 
subshells in the 660-740 eV range are compared to explore effects of the F ligands. The Xe 3d - ef continuum 
shape resonances dominate the photoionization cross sections of both the atom and molecule, but prominent 
resonances appear in the XeF 2 cross section due to nominal excitation of Xe 3d and F Is electrons to the lowest 
unoccupied molecular orbital (LUMO), a delocalized anti-bonding MO. The subshell ionization thresholds, the 
LUMO resonance energies and their oscillator strengths are calculated by relativistic coupled-cluster methods. 
Several charge states and fragment ions are produced from the atom and molecule due to alternative decay 
pathways from the inner-shell holes. Total and partial ion yields vary in response to the shape resonances and 
LUMO resonances. Previous calculations and measurements of atomic Xe 3d core-hole decay channels and 
our calculated results for XeF 2 guide interpretations of the molecular ion products. The partial ion yields of 
XeF 2 are dominated by Xe 3d core-hole decays, but distinct ion products are measured at the F Is - LUMO 
resonance. Xe 3d core-hole decays from XeF 2 produce lower charge states in comparison with atomic Xe, and 
energetic F ions are produced by Coulomb explosions of the molecular ions. The measurements support a 
model of molecular core-hole decay that begins with a localized hole, stepwise Auger electron emission spreads 
charge across neighboring atoms, and the system fragments energetically. 

PACS numbers: 32.80.Aa, 33.80.Eh 


I. INTRODUCTION 

The xenon fluorides, XeF„ (n = 2, 4, 6), have been 
the subject of numerous studies that have focused on 
their structures, physical, chemical, and spectroscopic 
properties. The nat ure of the chemical bonding in “no¬ 
ble gas compounds”fi^ or “compounds with an excess 
of electrons’’^! has been of continuing interest. Despite 
the closed-shell structure of atomic Xe, the high electron 
affinities of the F ligands induce bond formation by shar¬ 
ing charge with the Xe valence electrons. An early bond¬ 
ing model for linear XeF 2 suggested 3-center, 4-electron 
delocalized molecular orbitals (MOs) consisting of 2pz 
atomic orbitals (AOs) on the F atoms and the 5pz AO 
on Xe, where z is taken as the molecular axi^^. In this 
model, two electrons are in a bonding orbital delocal¬ 
ized over the three centers, and two electrons are in a 
nonbonding orbital on the F atoms only. The resulting 
picture is of a weakly-bonded system with mixed ionic 
and covalent character. The third, unoccupied 3-center 
orbital is the anti-bonding combination and was identi¬ 
fied as the lowest unoccupied molecular orbital (LUMO) 
in photoabsorption spectra.^^^ 

The Hartree-Fock self-consistent-held (SCF) calcula¬ 
tions of Basch et ol.® provide a more comprehensive 
model of bonding in the xenon fluorides and of their 
physical properties. Population analysis of the AO com¬ 
position of the MOs, the atomic charge distributions, 
and the orbital energies derived from the SCF calcu¬ 
lations are compared with electron spin resonance, nu¬ 


clear magnetic resonance, photoabsorption, and photo¬ 
electron measurements.® For example, core-level binding 
energy shifts measured in x-ray photoelectron spectra are 
sensitive to the valence electron charges on the Xe and 
F atoms.®^ More recent quantum chemical calculations 
treat relativistic and many-electron interactions in the 
xenon fluorides 

Here we explore molecular effects on inner-shell pho¬ 
toionization and core-hole decay by comparing measure¬ 
ments and calculations on XeF 2 with atomic Xe. Buth 
et explained how the observed increase in the decay 
width of Xe 4d holes with the number of F ligands results 
from interatomic electronic decay processes. Dunford et 
explored charge production, charge redistribution, 
and energetic ion fragmentation in the decay of Xe Is 
holes in XeF 2 . The initial core hole is highly localized, 
but as the core hole decays by the emission of fluores¬ 
cent photons and Auger electrons, holes are formed in 
the outer subshells, charge is redistributed to the F lig¬ 
ands, and the system undergoes Coulomb explosion. The 
kinetic energies of the F fragment ions measured by Dun- 
ford et al.ESI are smaller than the Coulomb energies that 
would be stored in the system at its ground state geom¬ 
etry. This suggests that ion fragmentation begins while 
the charge distribution is developing on the system. The 
picture that emerges is sketched in Fig. [^with core-hole 
decay, charge redistribution, and Coulomb explosion pro¬ 
ceeding concurrently on the femtosecond time scale. 

Inner-shell photoionization of molecules is also being 
explored with intense, ultrafast x-ray free-electron laser 
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FIG. 1: Illustration of inner-shell photoionization, core-hole 
decay, charge redistribution, and Coulomb explosion in 
XeF 2 . Inner-shell photoionization of Xe triggers core-hole 
decay that is initially localized but subsequently involves 
delocalized valence electrons. Charge spreads to neighboring 
atoms and the system Coulomb explodes. The processes 
proceed concurrently on the femtosecond time scale. 


(XFEL) pu lses t hat induce absorption of multiple x- 
ray photons. b^FS I recent pump-probe experiment on 
iodomethane, an optical laser dissociated the I atom and 
XFEL pulses triggered inner-shell photoionization and 
core-hole decay as a function of time delay.^ Charge 
transfer between the I and CH 3 fragments was measured 
as a function of internuclear separation, and the results 
are consistent with a classical over-the-barrier model of 
electron transfer. Such experiments provide new insights 
into molecular core-hole decay. 

In the present work, we used soft x rays in the 660- 
740 eV range to photoionize the Xe 3 d 5 / 2 , Xe 3 d 3 / 2 , 
and F Is subshells of Xe and XeF 2 . The Xe 3d and 
F Is ionization energies are similar, and the photoion¬ 
ization cross sections overlap, but the inner-shell states 
are atomic-like, so the initial core holes are localized. As 
Auger electron emission proceeds stepwise, delocalized 
valence electrons participate, charge is spread across the 
molecule, and the system fragments. An ion spectrome¬ 
ter was used to measure total and partial ion yields across 
the subshell edges. The 3d photoionization cross section 
of atomic Xe was measured previously,li^l^^ and varia¬ 
tions of the yields of the multicharged ions Xe'^+ (q = 
1-8) vs. photon energy were reported.^ The atomic Xe 
3d core-hole decay processes have been characterized us¬ 
ing Auger electron spectroscopy, calculated de-excitation 
pathways, and electron-ion and electron-electron coin¬ 
cidence spectroscopies q^j. gggj jg explore the 

effects of the F ligands on the photoionization cross 
sections and partial ion yields of XeF 2 . The Xe 3d - 


ef continuum shape resonances dominate the 3 d 5/2 and 
3 d 3/2 photoionization cross sections of both the atom and 
molecule, but distinct resonances are also observed in 
XeF 2 due to excitation of the LUMO from the Xe 3 d 5 / 2 , 
Xe 3 d 3 / 2 , and F Is subshells. The LUMO is a delocal¬ 
ized anti-bonding MO that was identified in early far- 
UV photoabsorption measurement^SHU ggd contributed 
to understanding the electronic structure of XeF 2 at the 
qualitative molecular orbital level. Theoretical treat¬ 
ments of inner-shell photoionization and core-hole decays 
in molecules are at an early stage compared with atoms. 
However, we performed relativistic coupled-cluster calcu¬ 
lations to determine the subshell ionization energies and 
the energies and oscillator strengths of the LUMO res¬ 
onances. The calculations help guide interpretations of 
the ion yields. 

Section |lT] of this paper describes the ion yield and 
branching ratio measurements, and computational meth¬ 
ods are discussed in Section IIIII In Section ||V] we dis¬ 
cuss the measured and calculated results for ion yields, 
branching ratios, energies and oscillator strengths. Con¬ 
clusions and suggestions for future research are in Sec. 

in 


II. EXPERIMENTAL METHODS 

Total and partial ion yields of Xe and XeF 2 were 
measured following photoionization of the Xe 3 d 5 / 2 , Xe 
3 d 3 / 2 , and F Is subshells in the 660-740 eV range. An 
ion time -of-flight spectrometer with pulsed extraction 
fielcPSM! was used on the variable line spacing plane 
grating monochromator (VLS-PGM) beamline at Wis¬ 
consin’s Synchrotron Radiation Center.l^ Xe was taken 
from a compressed gas cylinder, and the 5 mbar vapor 
pressure of XeF 2 was used from the solid at room temper¬ 
ature. Diffuse beams of the samples were produced with 
a capillary nozzle in the vacuum chamber and crossed 
with the photon beam in the interaction volume of the 
ion spectrometer. The ion spectrometer uses a pulsed 
electric field to push ions from the interaction region 
through an aperture into the time-of-flight drift tube for 
mass/charge dispersion and single-ion detection with a 
microchannel plate detector. The ion time-of-flight spec¬ 
tra from Xe and XeF 2 recorded over the 695-730 eV 
range are summed and plotted in Fig. The ion spectra 
in this photon energy range are dominated by Xe 3d core 
hole decays. The structures in the ion peaks are from 
the different masses and natural abundances of the Xe 
isotopes. The isotopically averaged mass of Xe is 131.3 
while F has a single isotope with mass 19.0. 

The ion spectrum of XeF 2 is plotted on an expanded 
scale in Fig. [^to show the weaker peaks. Strong Xe'^'*' 
(q = 1-6) peaks and weak XeF^, XeF^’’’, XeF+, and 
XeF^+ peaks were recorded, but the F+ (overlapped with 
Xe^"*") and F^’*' peaks are small. This is attributed to the 
low detection efficiency of the spectrometer for energetic 
ions produced by Coulomb explosions. The spectrometer 
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FIG. 2: Ion time-of-flight spectra of Xe (top) and XeF 2 
(bottom) summed over the 695-730 eV photon energy range. 

The structures in the peaks are due to the different Xe 
isotope masses. The weaker peaks of the XeF 2 spectrum are 
shown on an expanded scale in Fig. [3| 


is designed to measure the charge-state yields of atoms 
and par ent m olecular ions that have only small thermal 
energies.l^^^ Energetic fragment ions leave the interac¬ 
tion volume during delays between extraction pulses so 
most are not detected. The collection solid angles for en¬ 
ergetic ions are limited to a small range along the sym¬ 
metry axis of the spectrometer. The F+ and F^+ peaks 
measured over 695-730 eV are plotted in Fi^ The F+ 
and Xe^^ peaks overlap as shown in Fig. ®a) by com¬ 
paring with the Xe^"*" spectrum from atomic Xe. The F"*' 
ions are spread to higher and lower flight times than the 



time-of-flight (^s) 


FIG. 3: Ion time-of-flight spectrum of XeF 2 on an expanded 
scale to show the weaker peaks. The F+ and Xe^^ peaks 
overlap. The ion spectrometer has small detection 
efficiencies for energetic F+ and F^^ fragment ions. The F"*" 
and F^"*" spectra are also plotted in Fig. 


Xe^+ ions due to their momenta from fragmentation of 
molecular ions. The average kinetic energy of F+ is esti¬ 
mated to be ~15 eV from the peak splitting and extrac¬ 
tion field of the spectrometer. The splitting of the F^+ 
peaks in Fig. |^b) gives an average kinetic energy of ~27 
eV. Although energetic ions are not efficiently detected 
with our spectrometer, variations of the measured yields 
of F+ and F^+ give useful information and are discussed 
in Sec. IIVI 

Figure shows that the Xe charge state distribution 
shifts down in the molecule, which we attribute to redis¬ 
tribution of charge to the F ligands. For atomic Xe, the 
average charge over the 695-730 eV range is 4.5, while the 
average Xe charge for XeF 2 is 2.6 (see Table [lV| . We infer 
that, on average, two charges are redistributed to the F 
ligands. The parent ions, XeF^ and XeF^’’’, and the ions 
missing one F ligand, XeF+ and XeF^+, are weak com¬ 
pared with Xe'^+ (q = 1-6). This suggests that most of 
the molecular ions dissociate to three atomic ions. Since 
XeF 2 is a linear molecule and the two F ligands have the 
same mass, both much lighter than Xe, momentum con¬ 
servation shows that the Xe ions are produced with near 
zero velocities. Also, the Xe ion spectra show distinct 
isotope structure without broadening that would result 
from kinetic energy released in the fragmentation process. 
The Xe fragment ions are therefore efficiently extracted 
and detected, and the Xe'^’*' charge state distributions 
from the atom and molecule can be compared. 

For inner-shell photoionization, total ion yields vs. 
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FIG. 4: Ion spectra of and from photoionization of 
XeF 2 summed over the 695-730 eV photon energy range, 
(a) F^ overlapped with Xe^^ (blue curve) compared with 
Xe^"*" from atomic Xe (red curve). The F"*' component is 
split due to ~15 eV kinetic energy, (b) F^"*" peak splitting 
corresponds to ~27 eV kinetic energy. 


photon energy are proportional to photoabsorption cross 
sections to a good approximation and the fluorescence 
yield is small in our case. To measure total ion yields in 
this experiment, the pulsed extraction field was replaced 
with a static field that pushed all of the ions, includ¬ 
ing energetic fragment ions, into the spectrometer. The 
monochromator was calibrated by recording ion yields 
of Ne across the Is - 3p resonance at 867.12 eV^ and 
O 2 across the Is - resonance at 530.8 eV.I^ The 
monochromator resolution was ~0.7 eV at 700 eV. Ion 


spectra were recorded over ^660-740 eV with steps of 0.5 
or 1 eV. For each photon energy, the counts under each 
ion peak were summed and divided by the total counts 
to determine branching ratios. The total and partial ion 


yields vs. photon energy are discussed in Sec. IV 


III. THEORY 

In this section we present the method that we used 
to calculate the inner-shell excitations and binding en¬ 
ergies in XeF 2 at the Xe 3d and F Is edges. In our 
particular case, the goal was to calculate the transition 
energy of Xe 3d5/2 - LUMO, Xe 3d3/2 - LUMO, and 
F Is - LUMO, together with the corresponding bind¬ 
ing energy for Xe 3d5/2, Xe 3d5/2, and F Is. Molecu¬ 
lar inner-shell calculations are demanding; they require 
a high-degree treatment of both electron-correlation and 
orbital-relaxation effects, as strong relaxation of the wave 
function is induced by the removal of an inner-shell elec¬ 
tron. In addition, relativistic effects are more significant 
for inner-shell electrons than for valence electrons, since 
the former are in general moving faster. Here the explicit 
calculation of inner-shell electrons necessitates the use of 
all-electron relativistic methods. Therefore, in our calcu¬ 
lations we have employed all-electron relativistic coupled 
cluster (CC) approaches in order to obtain a systematic 
treatment aiming at high accuracy. 

ccP is a powerful numerical method used in quan¬ 
tum chemistry that utilizes an exponential parametriza- 
tion of the wave function in order to provide a rapidly 
convergent and size-extensive description for energy and 
properties. CC provides very good results for single¬ 
reference systems, whose electronic wave functions are 
dominated by a leading electronic configuration. In the 
CC framework, the calculations for excited or ionized 
states are based on the equation-of-motion (EOM) CC 
methods (also known as CC linear response theory in 
the literature)EOMCC is a “direct” approach that 
computes the difference between a target state (excited 
or ionized state) and a reference state (usually chosen as 
the closed-shell ground state), and thus can effectively 
exploit the cancelation of errors in the description of the 
two states. EOMCC calculations for core excitations or 
ionizations have been reported in the literature.l^SHHI As it 
is the spin-orbit coupling that introduces high computa¬ 
tional overheads for the fully relativistic CC calculations, 
in our present calculations we have adopted the cost- 
effective scheme of combining scalar-relativistic EOMCC 
calculations with spin-orbit corrections obtained at the 
Hartree-Fock (HF) level. Specifically, we have employed 
the spin-free exact two-co mpone nt theory in its one- 
electron variant (SFX 2 C-le)F 2 H 2 in combination with 
EOMCC methods, in order to achieve an efficient and 
accurate treatment of scalar-relativistic effects as well 
as a systematic incorporation of electron-correlation and 
orbital-relaxation effects. The SFX2C-le EOMCC re¬ 
sults for the excitation and binding energies are then aug- 
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merited with the spin-orbit contributions obtained at the 
level of Koopmans’ theorem using the Dirac-Coulomb- 
Gaunt (DCG) Hamiltonian ESI The non-relativistic HF, 
SFX2G-le HF, and SFX2C-le EOMCC calculations have 
been performed using CFOUR,E2lEIIwhile the DC and 
DCG HF calculations have been carried out using the 
DIRAG program package.E^HH 

In Table |l] we show the results for the binding energies 
calculated for Xe atom, HF molecule, and XeF 2 molecule. 
The equilibrium geometry calculated at the SFX2G- 
le/CCSD(T)/unc-ANO-RCC level (ixe-F = 1-9736 
has been used for the XeF 2 molecule, while the experi¬ 
mental bond length (rH-F=0.9168 A^has been used for 
the HF molecule. First, we present results obtained by 
Koopmans’ theorem using the non-relativistic, SFX2G- 
le, DC, and DCG Hamiltonians. The ANO-RCC basis 
in its uncontracted form (unc-ANO-RCCj^Tllll has been 
used in these calculations. The importance of relativis¬ 
tic effects is demonstrated by the difference between the 
non-relativistic and SFX2C-le results, while the spin- 
orbit corrections are obtained as the difference between 
the SFX2C-le results and those from DG and DCG cal¬ 
culations. In order to account for electron-correlation 
and orbital-relaxation effects, we perform several cal¬ 
culations at the SFX2G-le-EOM-GGSD and SFX2C-le- 
EOM-CCSDTEH levels. Systematically enlarged atomic- 
natural-orbital basis sets of double-zeta, triple-zeta, and 
quadruple-zeta quality (ANOO, ANOl, and AN02j^ 
as well as the unc-ANO-RCG basis have been used to 
demonstrate the basis-set convergence. In all the GG 
calculations, the Xe Is, 2s, and 2p orbitals have been 
frozen. Here it is interesting to observe the atomic na¬ 
ture for the inner-shell molecular orbitals: the correction 
energies for the Xe 3d (or the F Is) orbitals are similar 
in the Xe (or the F) atom than in the Xe (or the F) site 
of the XeF 2 molecule. The total binding energy has been 


obtained by combining the results of the best quality of 
accounting for various effects, including the DGG/unc- 
ANO-RCC results at the level of Koopmans’ theorem, the 
electron-correlation and orbital-relaxation corrections at 
the EOM-CCSD/unc-ANO-RGG level, and the triples 
corrections taken as the difference between the EOM- 
CCSDT/ANOl and EOM-CCSD/ANOl results. 

The positions of the LUMO resonance have been ob¬ 
tained by an additive scheme, by augmenting the accu¬ 
rately determined binding energies with the difference 
between the excitation and binding energies calculated 
at the SFX2G-le-E0M-CCSD/AN01 level (see the last 
row of Table |T| . In Table |h] we also show the calculated 
oscillator strengths for the corresponding resonant transi¬ 
tions in XeF 2 . We mention that, due to the D 2 h symme¬ 
try of the XeF 2 molecule as well as the atomic nature of 
the inner-shell orbitals, the atomic Xe 3d orbitals almost 
exactly map onto the corresponding molecular orbitals. 
Therefore we have used the atomic notation for the Xe 
3d orbitals infill The calculated and measured results are 


compared in Sec. IV 


Gonvergence problems in the solution of the EOMGG 
eigenvalue equations may occur for the standard David¬ 
son algorithm, when there exist multiply excited or 
“shake-up” configurations that are energetically quasi¬ 
degenerate with the target inner-shell excited or ion¬ 
ized states. Among our calculations, convergence diffi¬ 
culties have been encountered in the EOM-GGSD/unc- 
ANO-RCC and all the EOM-CCSDT calculations of the 
binding energies as well as the EOM-CCSD/ANOl cal¬ 
culations of excitation energies. Here we have adopted 
the Arnoldi algorithnJ^ to converge the solutions for 
these calculations. In the EOM-CCSD/unc-ANO-RGG 
calculations, virtual orbitals with orbital energies higher 
than 100 Hartree have also been excluded to expedite the 
convergence 



Xe 3d 

in Xe atom 

F Is in HF 

XeF 

2 molecule 


3d3/2 

3d5/2 

F Is 

Xe 3d3/2 

Xe 3d5/2 

F Is 

L 

koopmans’ theorem 

Nonrel/unc-ANO-RCC 

710.7 

710.7 

715.5 

714.4 

714.4 

715.0 

SFX2C-le/unc-ANO-RCC 

700.2 

700.2 

716.2 

703.9 

703.9 

715.6 

DC/unc-ANO-RCC 

708.1 

694.9 

716.3 

711.9 

698.6 

715.7 

DCG/unc-ANO-RCC 

707.5 

694.6 

716.1 

711.2 

698.3 

715.4 

SFX2C-le-EOM-CCSD corrections 

ANOO 

-10.9 

-10.9 

-19.2 

-11.5 

-11.5 

-20.2 

ANOl 

-16.3 

-16.3 

-20.0 

-16.7 

-16.7 

-21.0 

AN02 

-17.6 

-17.6 

-20.3 

-17.9 

-17.9 

-21.3 

unc-ANO-RCC 

-17.9 

-17.9 

-19.9 

-18.0 

-18.0 

-21.2 

EOM-CCSDT corrections on top of EOM-CCSD 

ANOO 

-0.9 

-0.9 

-1.7 

-1.2 

-1.2 

-2.9 

ANOl 

-0.7 

-0.7 

-1.8 

-1.0 

-1.0 

-2.9 

Total Binding Energy 

688.9 

676.0 

694.4 

692.2 

679.3 

691.3 

LUMO resonance 




682.6 

669.7 

683.1 


TABLE I: Computed binding energies and LUMO resonance positions (in eV). 
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F Is - LUMO 


Xe 3d 

- LUMO 



ls(F#l)+ls(F#2) 

ls(F#l)-ls(F#2) 

3d(2?5?VT 

3d{xz) 

3d(yz) 

3d(xy) 

3d(x^-y^) 

ANOO 

0.0627 

0 

0.0048 

0.0013 

0.0013 

0 

0 

ANOl 

0.0611 

0 

0.0047 

0.0013 

0.0013 

0 

0 


TABLE II: Oscillator strength for the F Is to LUMO and the Xe 3d to LUMO transitions in the XcF 2 molecule 

computed at the EOM-CCSD level. 


IV. RESULTS 
A. Total cross sections 


The total ion yield scans of Xe and XeF 2 are plotted in 
Fig.§ Our measurements of the Xe 3d cross sections of 

atomic Xe agree well with previous measurements 

but the XeF 2 cross section has not been reported pre¬ 
viously. Carroll et measured the binding energies 
of the Xe 3d5/2, Xe 3d3/2 and F Is subshells of Xe and 
XeF 2 , and our calculated binding energies match their 
measurements quite well (see Table III I. The nearly per¬ 
fect agreement between the computed and experimental 
numbers can partially be attributed to cancellation of 
remaining errors in the calculation, as the residual basis 
set error for the EOM-CCSD correlation correction and 
that for the triples correction appear to have different 
sign following the trend seen in Table The Xe dds /2 - 
3 d3/2 spin-orbit splitting is 12.8 eV in both Xe and XeF 2 , 
but the binding energies are shifted up by 2.9 eV in the 
molecule. The F Is binding energy of XeF 2 is shifted 
down by 5.5 eV relative to that of F 2 , which places it 
only 0.9 eV below the Xe 3d3/2 binding energy. The 
core-level binding energy shifts in XeF 2 are sensitive to 
the valence electron charge distributions on the Xe and 
F sites .1^ 


The atomic Xe 3d5/2 and 3d3/2 cross sections are dom¬ 
inated by strong, broad maxima near 688 eV and 702 
eV that are attributed to shape resonances in the f-wave 
continuum channels. The shape resonances can be un¬ 
derstood as independent-electron effects, (that is, those 
that can be understood in the orbital model of electronic 
structure) but measurement^^ and calculation^SMsH 
show coupling between the 3d5/2 - cf and 3d3/2 - ef chan¬ 
nels. Small features due to 6p and 7p Rydberg states 
were also recorded in the atomic X e cross section in agree¬ 
ment with previous measurements.1221221 

The Xe 3d - ef shape resonances produce strong max¬ 
ima near 688 eV and 700 eV in the XeF 2 cross sections, 
but distinct resonances also appear due to excitation of 
the LUMO. The Xe 3d5/2 - LUMO resonance is at 669.9 
± 0.3 eV, which is 9.4 eV below the 3d5/2 binding en¬ 
ergy and is in nearly exact agreement with theory. The 
calculated energies of the Xe 3d3/2 - LUMO and F Is - 
LUMO resonances are within 0.3 eV of each other, which 
is less than our 0.7 eV resolution, and we recorded a single 
peak at 682.8 ± 0.3 eV. However, the calculated oscillator 
strength of the F Is - LUMO resonance is ~ 20x larger 
than that of the Xe 3d3/2 - LUMO resonance (see Table 
[h]) Excitation of the LUMO was observed previously 


TABLE III: Binding energies and excited state energies 
from theory and experiment. 


target 

transition 

theory“ (eV) 

experiment (eV) 

Xe 

365/2 ^ 

676.0 

676.44 ± 0.05“ 

Xe 

363/2 ^ 

688.9 

689.23 ± 0.05“ 

HF 

F Is-i 

694.4 

694.25 ± 0.08*’ 

XeFa 

365/2 ^ 

679.3 

679.31 ± 0.05“ 

XeFa 

363/2 ^ 

692.2 

692.09 ± 0.05“ 

XeFa 

F Is"^ 

691.3 

691.23 ± 0.05“ 

XeFa 

3d5/2 “ LUMO 

669.7 

669.9 ± 0.3“ 

XeFa 

363/2 - LUMO 

682.6 

682.8 ± 0.3“ 

XeFa 

F Is - LUMO 

683.1 

682.8 ± 0.3“ 


“Ref. 1^ *'Ref. 1631 “This work. 


in far-UV photoabsorption spectra,in photoexcitation 
of the Xe 4d subshells,and in the Xe A-shell.tf^ 


B. Xe ion yields 


Ion spectra of atomic Xe were measured over 660-740 
eV and the Xe'^+ (q = 1-8) branching ratios are plotted 
in Figs. and The branching ratios strongly change 
as the photon energy increases from below to above the 
3 d5/2 ionization energy, particularly in response to the 
ef shape resonance. The relative yields of the lower 
charge states, q = 1-3, decrease while q = 4-7 increase. 
Higher charge states are expected, as photon energies ac¬ 
cess deeper subshells and open additional core-hole decay 
steps.l^ Smaller changes in the branching ratios appear 
as the photon energy increases above the 3d3/2 ionization 
energy and excites its shape resonance. A second maxi¬ 
mum also appears in the 3d5/2 cross section at ~707 eV 
due to coupling with the 3d3/2 shape resonance 

The present measurements of Xe'^“'" branchi^ ratios 
are compared with those of Saito and SuzukflS in Ta¬ 


ble IV for two energy regions. The 660-670 eV region 
is below the 3d5/2 - 6p Rydberg resonance at 674 eV, 
which is the lowest excitation energy of a 3d electron. 
The ion branching ratios in that region are determined 
by photoionization of the 5s and 5p valence electrons and 
populate lower charge states compared with photon en¬ 
ergies that excite 3d electrons. Small variations of the 
branching ratios are observed in Figs. and from ex¬ 
citation of the 3d5/2 - 6p Rydberg state at 674 eV, but 
the ion yields are clearly dominated by the 3d - ef shape 
resonances. 


Table IV also compares measurements made in the 
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TABLE IV: Branching ratios for Xe ion charge states 
from Saito and SuzukP^ at 670 eV and 700 eV 
compared with the present results for Xe averaged over 
660-670 eV and 695-730 eV. The last column shows the 
present results for Xe ion charge states from XeF 2 
averaged over 695-730 eV. 


ion 670 eV“ 660-670 eV*’ 700 eV“ 695-730 eV** XeFz: 695-730 eV'’ 


Xe+ 

3.5 

4.9 ± 0.2 

0.7 

0.7 ± 0.1 

19.4 ± 

Xe^"^ 

29 

32.4 ± 0.3 

4.0 

5.5 ± 0.1 

35.7 ± 

Xe®"^ 

25 

23.3 ± 0.2 

7.5 

7.8 ± 0.1 

25.6 ± 

Xe4+ 

30 

26.4 ± 0.3 

41 

39.8 ± 0.2 

12.1 ± 

Xe®"^ 

9.0 

7.8 ± 0.1 

29 

27.9 ± 0.2 

4.2 ± 

Xe«+ 

3.0 

3.2 ± 0.1 

15 

15.0 ± 0.1 

1.7 ± 

Xe*”"^ 

0.5 

1.3 ± 0.1 

2.5 

2.8 ± 0.1 

1.3 ± 

Xe®"^ 

0.1 

0.6 ± 0.1 

0.3 

0.5 ± 0.1 

- 

Xe®+ 

- 

0.1 ± 0.1 

- 

- 

- 



■"Ref. HU 

'’This work. 


695-730 

eV 

range above 

the 3d5/2 

and 3d3/2 

binding 


energies. Ionization of a 3d electron triggers a vacancy 
cascade, a series of Auger electron transitions through 
alternative pathways that lead to the observed range of 
charge states.l^SIlZlThe de-excitation pathways have been 
calculated using codes for the electronic s truct ures of 
the hole states and Auger transition rates.l^^*^ Corre¬ 
sponding calculations for core-hole decays in molecules 
are more chall engin g due to their non-spherical, multi¬ 
site structures.l^^Eni However, the atomic measurements 
and calculations provide guidance to understanding the 
first decay steps of Xe 3d holes in both Xe and XeF 2 . 


C. XeF2 ion yields 

The Xe 3 d 5/2 - ef and Xe 3 d 3/2 - ef shape resonances 
dominate the total ion yield scans of Xe and XeF 2 plotted 
in Fig. Hartree-Fock-Slater calculations of the Xe 
3d and F Is photoionization cross sections averaged over 
700-750 eV give a Xe 3d cross section 7x larger than 
twice the F Is cross section.l^ The only clear signature 
of F Is excitation is at the 683 eV LUMO resonance 
that also has a contribution from Xe 3 d 3/2 excitation 
(see Tables HB- Away from the 683 eV resonance, Xe 
3d core hole decay is the dominant contributor to the 
ion yields of XeF 2 . Previous work on atomic Xe suggests 
a model for XeF 2 in which a Xe 3d hole decays to two 
4s, 4p, or 4d holes that are localized on followed 

by Auger transitions involving delocalized Xe 5p -I- F 
2p MOs that spread charge to the F sites. The shift to 
lower Xe'^'*' charge states in going from the atom to the 
molecule shown in Fig. is quantified in Table |IV| by 
comparing measurements over 695-730 eV. With charge 
spread across the F and Xe sites, the system Coulomb 
explodes. Energetic fragmentation is apparent in the ion 


0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 


Jonauskas et reported that 90% of 3d holes decay 
to double-hole states in the N shell, i.e., 4d“^, 4p“^4d“^, 
4s“^4d“^, 4p“^, and 4s“^4p“^. Those holes are still lo¬ 
calized. Only 4% of 3d holes decay to states with 5s or 
5p holes, i.e., 4d“^5p“^ and 4d“^5s“^. It is not until 
the second Auger decay step that a significant fraction 
of 5p holes are produced.^ In XeF 2 , the Xe 5p electrons 
combine with F 2p electrons to form delocalized valence 
MOs.l^ By comparison with the stepwise 3d hole decays 
in atomic Xe, we infer that in XeF 2 delocalized electrons 
participate in the second Auger decay steps of a Xe 3d 
hole and spread charge to the F sites. 


spectra of F+ and F^+ plotted in Fig. We explained 
in Section [n] that our ion spectrometer has relatively low 
collection efficiencies for energetic F+ and F^+ ions but 
good sensitivity for Xe’*’*' ions. We can therefore compare 
Xe'*"'" ion yields from the atom and molecule. In addition, 
although the measured F"*" and F^"*" yields are small, their 
variations with photon energy are informative. 

Ion spectra of XeF 2 were recorded over 665-730 eV 
and the counts under each peak summed and divided 
by the total counts to determine branching ratios. No 
corrections were made for the low collection efficiencies 
of F"'" and F^"*". The variations with energy are plotted in 
Figs. [8}|T0| Features are observed in response to the Xe 
3 d 5/2 - LUMO resonance near 670 eV, at the unresolved 
Xe 3 d 3/2 - LUMO and F Is - LUMO resonances near 683 
eV, and from the two Xe 3d - ef shape resonances with 
maxima at 688 and 700 eV. 

The first excitation of a Xe dds /2 electron occurs at the 
670 eV LUMO resonance that is located 9.4 eV below the 
3 d 5/2 ionization energy. The calculations of BasclP char¬ 
acterize the LUMO as an anti-bonding MO with Mulliken 
populations of 68 % Xe 5p, 31% F 2p, and 1% F 2s. We 













expect the LUMO electron to act as a spectator in the 
first Auger decay step of the 3 d 5/2 hole to two holes in 
the Xe N shell but to participate in subsequent decay 
steps. Our branching ratio measurements show decreases 
of Xe+, Xe5+, Xe®+, XeF+, XeFl+, XeF+, and XeF^+. 
The branching ratios for Xe^+, Xe^"*", Xe"^'*', and F^+ in¬ 
crease, while the overlapped F+ and Xe^+ ions shows no 
response. 

The most dramatic variations of the branching ratios 
occur at the 683 eV resonance where the LUMO is ac¬ 
cessed by both F Is and Xe 3 d 3/2 electrons. The Xe"*", 
F+, and XeF+ ratios increase while the Xe®'*' (q = 2- 
6 ) and F^"*" ratios drop. The ion spectrum recorded on 
resonance is plotted in Fig. [TTIand shows enhanced pro¬ 
duction of Xe"*" and F+ compared with Fig. In com¬ 
parison with Fig. I^a), the F"*" ions are much stronger 
on resonance than the Xe^"*" component. This is con¬ 
sistent with the calculated oscillator strengths in Table 
|TT] that find the F Is - LUMO resonance to be ^ 20 x 
stronger than the Xe 3 d 3/2 - LUMO excitation.l^ The F 
Is - LUMO transition moment is parallel to the molec¬ 
ular axis, and the photon beam is highly linearly polar¬ 
ized (0.996 ± 0.003).!^ Resonant excitation selectively 
excites molecules aligned with the polarization direction. 
The ion spectrometer axis is aligned with the polariza¬ 
tion, which contributes to the measured enhancement of 
the F+ yield. The splitting of the F+ peak is smaller 
on resonance and corresponds to a fragmentation energy 
of only ^3 eV. These results distinguish F Is excitation 
and decay from Xe 3d excitation and decay that dom¬ 
inates at other energies. F Is photoionization of sev¬ 
eral small molecules (SFg, SiF 4 , BF 3 , CF 4 ) primarily 
produces singly- and doubly-charged fragment ions.^^®^ 
This is consistent with a simple picture of F Is photoion¬ 
ization followed by a single Auger decay step that pro¬ 
duces low total charge states. With low total charge on 
the molecular ion, the kinetic energies of the fragment 
ions are smaller than for decays of Xe 3d holes. 

The Xe 3d shape resonances dominate at photon en¬ 
ergies above the 683 eV resonance. The branching ra¬ 
tios of higher Xe'^+ charge states (q = 3-6) and F^+ 
increase while the ratios of lower charge states (Xe’*', 
XeF^, XeF 2 ~'', XeF+, and XeF^+) decrease. The F+ 
yield in this energy range is unclear due to overlapping 
with Xe^+. However, enhancement of the F^+ yield and 
the high kinetic energies of F+ and F^+ estimated from 
the ion spectra in Fig. |^are characteristic of Xe 3d hole 
decay in XeF 2 . 


V. CONCLUSION 

Inner-shell photoionization and core hole decay of Xe 
3 d 5 / 2 , Xe 3 d 3 / 2 , and F Is subshells of Xe and XeF 2 are 
compared over 660-740 eV to explore effects of the F 
ligands. The total and partial ion yields vary with pho¬ 
ton energy in response to Xe 3 d 5/2 - ef and Xe 3 d 3/2 - 
ef shape resonances in both the atom and molecule as 


well as excitation of the XeF 2 LUMO from Xe 3d and F 
Is subshells. Coupled-cluster calculations accounting for 
relativistic effects give accurate values for ionization en¬ 
ergies and resonance energies, and calculated oscillator 
strengths have helped to guide interpretations of mea¬ 
surements at the LUMO resonances. 

We use the stepwise Auger decay model for 3d holes in 
atomic X^^ to propose a corresponding model for XeF 2 . 
The first Auger decay step produces two holes in the Xe 
N shell that are localized on Xe, but the following Auger 
transitions spread charge across the three atomic sites. 
The system Coulomb explodes, producing Xe'^''" (q = 1- 
6 ) and energetic F+ and F^+ ions. The molecule responds 
differently to excitation of the F Is - LUMO resonance. 
The relative yields of F+ and Xe'*' increase, the F^+ and 
Xe'^+ (q = 2 - 6 ) yields decrease, and the fragmentation 
energy of F+ ions is relatively small. Therefore, XeF 2 
molecule is very exceptional in the photon-energy range 
studied. First, the molecular core-hole decay can be com¬ 
pared with Xe atom, and second, either an electron from 
the F site or the Xe site can be excited/ionized giving 
rise to a totally different core-hole decay. 

Additional measurements are suggested by the results. 
Momentum-resolved, multiple-ion-coincidence measure¬ 
ments with high-resolution synchrotron radiation can di¬ 
rectly c orrela te the F and Xe fragment ions and their 
energiesCore-hole decay and molecular fragmen¬ 
tation proceed on the femtosecond time scale, which is 
much faster than synchrotron x-ray pulses. However, the 
intense femtosecond x-ray pulses produced at XFELs can 
be used for time-resolved measurements.li^ New capabil¬ 
ities are being developed at XFELs, such as using two 
femtosecond x-ray pulses with different colors and vari¬ 
able time delay.^ The first pulse can trigger a core hole 
decay process in a molecule while the second pulse probes 
the dynamics . 1 ^ 

Core hole decay dynamics in molecules is a challenging 
problem for theory and computational methods. Treat¬ 
ments require a proper description of orbital relaxation 
and electron correlation. Moreover, inner-shell exci¬ 
tations induce multiply excited configurations due to 
“shake-up” processes making the convergence of excited 
state calculations very difficult. The lineshape of inner- 
shell excitations can be elucidated by calculating the 
Auger decay of the inner-excited states. Calculations of 
Auger decays need an accurate description of the contin¬ 
uum for the emitted Auger electron. This is quite chal¬ 
lenging to implement within commonly used bases, such 
as Gaussian bases, and some m ethod s have been imple¬ 
mented to address this problerrPUlll. In our particular 
case, the decay from Xe 3d holes in XeF 2 , two or more 
Auger transitions need to be treated while the molecu¬ 
lar geometry expands and fragments, involving also an 
important motion of the nuclei. 
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FIG. 5: Total ion yield scans of Xe and XeF2. Panel (a) is 
the Xe scan with Rydberg resonances marked by red lines: 
3d5/2 - 6p (673.9 eV), 3d3/2 - 6p (686.6 eV), and 3d3/2 - 7p 
(688.2 eV). The blue lines mark the 3d5/2 (676.44 eV) and 
3d3/2 (689.23 eV) binding energies.® The strong, broad 
peaks are the Xe 3d5/2 - rf (~688 eV) and Xe 3d3/2 - rf 
(~702 eV) shape resonances. Panel (b) is the XeF2 scan 
with resonances marked by red lines: Xe 3d5/2 - LUMO 
(669.9 eV) and the unresolved Xe 3d3/2 - LUMO and F Is - 
LUMO resonances (682.8 eV). The blue lines mark the Xe 
3d5/2 (679.31 eV), Xe 3d3/2 (692.09 eV), and F Is (691.23 
eV) binding energies.® The strong, broad peaks are the Xe 
3d5/2 - ef (~688 eV) and Xe 3d3/2 - d (~700 eV) shape 
resonances. 
















12 



I ^_ I _^^_u 

660 680 700 720 740 

energy (eV) 


FIG. 6: Branching ratios for (q = 1-4) from 

photoionization of Xe. The vertical lines mark the 3d5/2 
(676.44 eV) and 3d3/2 (689.23 eV) binding energies. 
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FIG. 7: Branching ratios for (q = 5-8) from 
photoionization of Xe. The vertical lines mark the 3d5/2 
(676.44 eV) and 3d3/2 (689.23 eV) binding energies. 
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FIG. 8: Branching ratios for Xe+, F+/Xe^+, and 

from photoionization of XeF2. The blue vertical lines mark 
the Xe 3d5/2 (679.31 eV), Xe 3d3/2 (692.09 eV), and F Is 
(691.23 eV) binding energies. The red vertical lines mark 
the Xe 3d5/2 - LUMO (669.9 eV) and Xe 3d3/2,F Is - 
LUMO (682.8 eV) resonances. 
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FIG. 9: Branching ratios for (q = 3-6) from 
photoionization of XeF2. The blue vertical lines mark the 
Xe 3d5/2 (679.31 eV), Xe 363/2 (692.09 eV), and F Is 
(691.23 eV) binding energies. The red vertical lines mark 
the Xe 365/2 - LUMO (669.9 eV) and Xe 363/2 ,F Is - 
LUMO (682.8 eV) resonances. 
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FIG. 10: Branching ratios for XeF’*", XeF^"*", XeFj and 
XeF2^ from photoionization of XeF2. The blue vertical lines 
mark the Xe 3d5/2 (679.31 eV), Xe 3d3/2 (692.09 eV), and F 
Is (691.23 eV) binding energies. The red vertical lines mark 
the Xe 3d5/2 - LUMO (669.9 eV) and Xe 3d3/2,F Is - 
LUMO (682.8 eV) resonances. 
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FIG. 11: Ion spectrum of XeF2 recorded on the Xe 3d3/2, F 
Is - LUMO (~683 eV) resonance. The relative yields of Xe"*" 
and F"*" fragment ions increase on resonance. The splitting 
of the F^ peak corresponds to ~3 eV kinetic energy. 




